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Abstract 

The quantum theory of optical coherence is apphed to the scrutiny of the statistical 
properties of the relic inflaton quanta. After adapting the description of the quantized scalar 
and tensor modes of the geometry to the analysis of intensity correlations, the normalized 
degrees of first-order and second-order coherence are computed in the concordance paradigm 
and are shown to encode faithfully the statistical properties of the initial quantum state. The 
strongly bunched curvature phonons are not only super-Poissonian but also super- chaotic. 
Testable inequalities are derived in the limit of large angular scales and can be physically 
interpreted in the light of the tenets of Hanbury Brown-Twiss interferometry. The quantum 
mechanical results are compared and contrasted with different situations including the one 
where intensity correlations are the result of a classical stochastic process. The survival 
of second-order correlations (not necessarily related to the purity of the initial quantum 
state) is addressed by defining a generalized ensemble where super-Poissonian statistics is 
an intrinsic property of the density matrix and turns out to be associated with finite volume 
effects which are expected to vanish in the thermodynamic limit. 
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1 Formulation of the problem 



The current data accounting for various properties of the Cosmic Microwave Background 
(CMB) anisotropics and polarization are all consistent with the presence, prior to matter- 
radiation equality, of a quasi-flat spectrum of curvature perturbations whose origin is cus- 
tomarily attributed to an (early) inflationary stage of expansion occurring when the Hubble 
rate was, roughly, of the order of one millionth of the Planck energy scale (see Refs. [H |2l |3l H] 
and, for a more comprehensive perspective on large-scale data, also Ref. 0). When CMB 
data are combined with other data sets (such as the large-scale structure data P, [7j and the 
supernova data P |9]) the typical parameters describing the large-scale curvature modes are 
slightly (but not crucially) modified. The agreement between the parameter determinations 
obtained by combining different data sets seems to reach the level of the per mill. This 
apparent (statistical) accuracy can be reached within the simplest scenario, conventionally 
dubbed ACDM where A stands for the dark energy component and CDM for the cold dark 
matter component. The quest for statistical accuracy should not forbid to think, in broader 
terms, to the very nature of the pre-infiationary initial conditions whose precise nature is all 
but well established [5] and anyway not explained in the framework of the ACDM paradigm. 

Can we establish, independently of the CMB data sets (see e.g. Refs. [U, Ej HI |3]), the 
duration of the inflationary stage of expansior§]? Not really: inflation could have been very 
long or it could have been just minimal (i.e. with approximate duration between 63 and 65 
e-folds by assuming the largest value of the slow-roll parameter compatible with the current 
upper limits set by the WMAP data [HE])- Is it known which were the initial conditions of 
the scalar and tensor modes of the geometry around the onset of the inflationary dynamics? 
Not really: if inflation lasted much more than the required 65 e-folds probably the only 
sound initial state for the inflaton quanta (and within the logic of inflationary models) was 
the vacuum. Conversely, if the duration of inflation was just minimal (or close to minimal) 
then different kinds of quantum mechanical initial states could play a decisive role and their 
associated energy density can only be constrained by back- reaction considerations |10j . 

The pair of questions contained in the previous paragraph are usually answered within 
two opposite points of view. Within the first set of hypotheses, the duration of the infla- 
tionary phase is determined by the aim of reproducing some structure in the temperature 

^In the ACDM paradigm, the total curvature of the Universe receives, at the present time, a leading 
contribution from the extrinsic curvature and a subleading contribution from the intrinsic (spatial) curvature. 
The role of inflation is, in this context, to make the ratio between intrinsic and extrinsic curvature sufficiently 
minute at the onset of the radiation dominated epoch so that it can easily be of order 1 today (during a 
decelerated stage of expansion the ratio between intrinsic and extrinsic curvature is actually increasing). 
The duration of the inflationary phase required to solve the latter goal represents one of the ways of pinning 
down the minimal number of inflationary efolds. 
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or polarization angular power spectra: in this case the duration of the inflationary phase is 
correlated (by construction) with the spectral behaviour of the initial data for the scalar and 
tensor modes of the geometry. The second (often tacit) assumption is that the duration of 
inflation is much larger than 60 efolds: in this way quantum mechanical initial conditions 
are better justified but the initial state does not have specific observable consequences. In 
different terms, the present measurements on the large-scale temperature and polarization 
anisotropies are only able to probe the power spectrum of curvature (adiabatic) fluctuations 
possibly present prior to matter-radiation equality over typical wavelengths much larger than 
the Hubble radius at the corresponding epoch. The increase of the number of inflatons or the 
addition of supplementary (non adiabatic) components in the initial conditions [121 lEl [E] 
automatically increases the number of parameters by making the model less predictive even 
if, potentially, more sound. 

The aim of the present paper is less pretentious than the two extreme approaches men- 
tioned in the previous paragraph: instead of arguing (on a purely theoretical basis) how long 
inflation must have been, it seems plausible to scrutinize wether it is possible, at least in 
principle, to determine the statistical properties of the initial state of relic inflaton quanta. 
In this respect it is both plausible and useful to draw a physical analogy with a similar class 
of problems arising in the quantum optical treatment of the fluctuations of visible light (see 
p3| \T6\ [T7] for three classic treatises covering, with different emphasis, all the theoretical 
tools which will also play a role in the forthcoming considerations). Consider, for sake of 
concreteness, a scalar quantum field V{x,t) where x denotes the spatial coordinate and r 
the time variable; V{x, r) might denote the quantum field describing either relic phonons 
or the single polarization of a graviton; in the following discussion, however, the scalar field 
V{x, t) denotes the single polarization of an electromagnetic wave in the visible frequencies 
as sometimes done in quantum optics [15]. The correlation functior§|: 

G'«(x,y;n,r2) = {V{x,r,)V{y,T^)) (1.1) 

can be probed, in quantum optics, by the Young two-slits experiment which is sensitive to the 
interference of the amplitudes of the radiation field [17] (see Fig. [T]). Young interferometry 
is not able, by itself, to provide information on the statistical properties of the quantum 
state of the radiation field since various states with diverse physical properties (such as laser 
light and chaotic light) lead to comparable degrees of first-order coherence. The normalized 

■^The averages (...) can denote either ensemble (statistical) or quantum averages depending on the nature 
of the source and also upon the preferred physical description. In the present paper the amplitude and the 
intensities are always related to quantum fields unless explicitly stated (see, in particular, section [6]). 
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counterpart of Eq. f ll.ip can be written as 



{x,y]Ti,T2) 



{V{x,n)V{y,T2)) 



(1.2) 
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^{\V{x,nW){\V{y,r,W)' 



defines the degree of first-order coherence. In Young interferometry the hght field can be 
classified depending upon the value of g^^\x, y; ti, T2): 

• the light is first-order coherent provided g^^\x,y]Ti,T2) = 1; 

• the light is partially coherent is < g^^\x, y; ri, T2) < 1; 

• the light is incoherent if g^^\x, y; ti, T2) = 0. 

In the Young two-slit experiment the electric fields emerging from the two pinholes produce 
the interference fringes on the second screen. The maximal (total) intensity on the second 
screen can be written as Ima.x = Xi +X2 + 2^/T^T^g^^\T) while the minimal (total) intensity 
on the second screen can be written as Xmin = Xi + X2 — 2^yZlX2g^^\T) where Xi and X2 
denote the intensity of the radiation field in each of the two pinholes. The visibility, i.e. 
(Xmax — Xmin)/(Xmax + Xmin) coiucidcs cxactly with g^^^r) in the case Xi = X2. If^?«(r) = l 
the visibility is maximized and, as discussed before, the light is said to be first-order coherent. 

Until the mid fifties, Eq. (11.11) has been used to define the coherence of the radiation field: 
a field was said to be coherent when the interference fringes are maximized in Young-type 
(two-slit) correlation experiment of the type of the one reported in the left plot of Fig. [H 
The applications of the Hanbury Brown-Twiss (HBT) effect (first to stellar interferometry 
[T8] and then more specifically to quantum optics [19]) demanded the accurate study of not 
only the correlations between field strengths (as defined in Eq. (11.11) ) but also the analysis 
of the correlations between the intensities of the radiation field, i.e. 



where now X(x, r) = V'^{x, r) is the intensity of the radiation field. In the realistic case when 
the radiation field is described by the electric field the intensity will simply be the squared 
modulus of the electric field. The physical implications of HBT interferometry in general 
and of Eq. (II. 3p in particular have been important for many areas of physics ranging from 
stellar astronomy [18] and quantum optics [IHl [HI |17], to pion interferometry [201 EU [22] 
and subatomic physics (see [231 El] for two comprehensive reviews). In subatomic physics 
HBT interferometry has been used to determine the hadron fireball dimensions [22] which 
is related to the linear size of the interaction region in proton-proton collisions. In Fig. [1] 
Young interferometry (plot on the left) is schematically compared to HBT interferometry 



G^'^\x,y;Ti,T2) = {i{x,Ti)i{y,T2)) 



(1.3) 
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Figure 1: The Young (plot on the left) and Hanbury Brown-Twiss (plot on the right) exper- 
iments are compared; r denotes the time delay (the speed of light c has been restored while 
natural units h = c = 1 will be used throughout). 



(plot at the right). In short the difference between the two experiments resides in the order 
of the correlation. In the case of Young interferometry (plot on the left in Fig. [1]) the 
fringes on the second screen arise as the result of the interference electric field amplitudes. 
In the case of HBT interferometry the intensities of the radiation field are measured at the 
correlator (see right plot in Fig. [T]). Since the intensities are quadratic in the amplitude of 
the electric fields, HBT involves the study of second-order correlation effects. Furthermore, 
since in quantum theory the intensities of the radiation field are quantized, HBT correlations 
represented some of the first evidence of quantum effects in the description of optical fields 
[15] as neatly expressed by Glauber [25], [26| (see, in particular, section IV of the first paper 
quoted in Ref. [26]). 

The Glauber theory of optical coherence [251 126] generalizes the concept of first-order 
coherence to higher orders and, in particular, to second-order. The second-order correlator 
defined in Eq. fll.3p can be written in its normalized form in full analogy with Eq. fll.21) . In 
the case when the intensities are purely classical stochastic variables the degree of second- 
order coherence can be written as 



g^'^\x,y;Ti,T2) 



{X{x,Ti)Z{y,T2)) 



:i.4) 



It is also common to define a slightly different normalized correlator (see, e.g. [23] 

(X(f,r)X(y,r)) 



R{x, y, r) 



1. 



:i.5) 



(X(f,r))(X(y,r)) 

If the intensities are constructed from an appropriate field operator, Eq. f ll.4p is usually 
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written, in a quantum optical context, as 

\x,y-n,r2) = . . (1.6) 

(: X(x,ri) :) (: X(y,r2) :) 

where the colon makes explicit the normal ordering of the operators. In quantum optics 
it is natural to impose the normal ordering in the correlators since the detection of light 
quanta (i.e. in the optical range of frequencies) occurs by detecting a current induced by 
the absorption of a photon [251 ES]. In Fig. [T] (plot at the right) the basic logic of the HBT 
experiment is schematically illustrated: the electric field is first split into two components 
through the beam splitter, then it is time-delayed and finally recombined at the correlator. 
The HBT setup provides therefore an operational definition for correlating the intensities of 
the radiation field. Conversely Young interferometry only probes the correlations between 
the amplitudes of the radiation field. According to the Glauber theory of optical coherence, 
the radiation field is said to be first-order coherent if g^^\x, y; ri, T2) = 1; the radiation field 
is said to be second-order coherent if g^'^\x, y; ti, T2) = 1. As we shall see in a moment, the 
coherent states of the radiation field are both first-order and second-order coherent. 

HBT interferometry encodes two complementary pieces of information characterizing the 
source, i.e. 

• the linear (or angular) size of the emitting (hyper) surface; 

• the statistical properties of the emitting quanta (photons, pions, phonons, gravitons). 

The first aspect is illustrated in the left plot of Fig. [2] where the correlations of the intensities 
of the radiation field X(x, ri) and X{y,T2) are depicted in the situation where ti = T2 = t. 
In Eqs. (II. 4p and (11. 5p {Z{x,T)X{y,T)) denotes the intensities measured both in x and in 
y while (X(x, r)) and (X(y, r)) denote the intensities measured separately in the two points. 
This definition is schematically illustrated in the left plot of Fig. [2l In a quantum mechanical 
perspective Eqs. (I1.5p - (I1.6I) imply that the normalized degree of second-order coherence can 
be measured by counting the photons. More specifically, supposing that ri = T2 the number 
of particles observed simultaneously in x and y can be divided by the product of the number 
of counts observed separately in x and y (see also Fig. [2]). The same reasoning holds, 
of course, also in the case of pions [22] as well as in the case of other particles obeying 
the Bose-Einstein statistics such as relic phonons and relic gravitons. In this sense HBT 
interferometry is deeply connected to what is known, in high-energy physics, as the study of 
Bose-Einstein correlations [22l [23] . 

In the case of the relic phonons and gravitons the normal ordering of the correlators is 
not specifically justified (even if it is technically useful, as we shall see). One of the purposes 
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Figure 2: In the plot on the left the physical idea of interfering intensities is schematically 
illustrated. In the plot at the right the different values of the intercept ng'^{0) are reported 
for different quantum states as a function of the average multiplicity of each quantum state. 

of the present study is to give the correct quantum mechanical definition of the degree of 
second-order coherence in the case of the relic phonons and of the relic gravitons in strongly 
correlated quantum states with large occupation numbers per Fourier mode. 

If the detection of the intensity occurs at the same spatial location (i.e. x = y) the degree 
of second-order coherence will depend upon the time difference t = ti — T2 with a bell-like 
shape. For r — )■ (zero time-delay limit) the degree of second-order coherence reaches a 
specific value which depends upon the statistical properties of the source and which will be 
discussed in a moment. When g^'^^r) — 1 in the limit r ^ 1, the width of g{T) estimates 
the coherence time of the source. A similar discussion can be conducted for the degree 
of space-time coherence and allows, in pion physics, an approximate determination of the 
hadronic fireball dimensions [22] (see also [231 121] ) • 

If photons are detected by photoelectric counting [TB], Eq. fll.6p can be written, for a 
single of the radiation field and in the limit r — (i.e. zero time-delay limit) as 

5C^.(0) = = + 1. D'^ = (AT^) - (iV)^ (1.7) 

where a and obey the usual Heisenberg-Weyl algebra [a, d'''] = 1 and D'^ denotes the 
variance. To pass from the first equality in Eq. (11. 7p to the second expression involving 
it must be noted that, using the commutation relations, (a'^a^aa) = (iV^) — (N) where 
N = d^a. Using now the definition of we also have that (d'^d^ad) = [D"^ — (N) + (iV)^]. 
The second equality in Eq. (I1.7P is finally proven by appreciating that, in the definition of 
^*^^^(0), the term {a^ aa) is divided by (a^a)'^. 



7 



From Eq. (11. 7p it is clear that different quantum states will lead to different values 
of g^^\0). For practical applications it is useful to define the so-called Mandel parameter 
(related to the zero-time-delayed correlator of Eq. (11.51) ) whose expression is 

Q=(iV)[^(2)(0)-l] = ^-l. (1.8) 



In the case of a coherent state 

D^ = {N), g^^\0) = l, and Q = 0; (1.9) 

and the radiation field is said to be second-order coherent. If the quantum state coincides 
with a Fock state containing n particles (i.e. |n)), the Mandel parameter equals —1. The 
(single mode) Fock states defines the lower limit of the degree of second-order coherence 
i.e. Q > —1 and ^'•^^(0) > 1 — 1/ {N) (the equality is reached exactly in the case of a Fock 
state). Chaotic (e.g. white) light has the property of leading to a degree of second-order 
coherence double than in the case of a coherent state, i.e. ^'•^^(0) = 2 which is a direct 
consequence of the (single mode) density matrix for a thermal state (see e.g. [IS]). The 
quantum mechanical correlations are then reflected in the degree of second-order coherence 
and, ultimately, in the magnitude and sign of ^^^^(0). In Fig. [2] (iV)^^(O) is illustrated for 
different quantum states as a function of the average multiplicity n = (N) of each state. To 
distinguish graphically the different states it is practical to plot ng^'^\0) as a function of n 
(rather than 'g^'^^ (0) itself as a function of n) . Chaotic light is an example of bunched quantum 
state (i.e. ^'^^■'(O) > 1 implying more degree of second-order coherence than in the case of a 
coherent state). Fock states are instead antibunched (i.e. ^*^^-*(0) < 1) implying a degree of 
second-order coherence smaller than in the case of a coherent state. Experimentally the zero 
time-delay limit is justified because the counting of photons (or pions) is made for typical 
times smaller than the coherence time of the source. In this sense bunched particles tend to 
arrive at the photodetector more simultaneously than their antibunched counterpart. The 
concept of bunching will be relevant for a complete understanding of the physical properties 
of curvature phonons (see, e. g. sections [2] and E]) . 

In this paper it will be argued that the tenets of the quantum theory of optical coherence 
can be used to fully characterize the correlation properties of cosmological perturbations. 
To explore different sets of initial conditions in conventional inflationary models the idea has 
been often to play either with a pre-infiationary phase of limited duration or to assign, in 
fully equivalent terms, an initial state on a given space-like hypersurface (see e.g. 1271128] 
and references therein). Different initial states result in large-scale modifications of the 
power spectrum which can be used either to suppress or to increase the power at large 
scales [28] (see also [29l [30]). Still, as it will be shown, different initial state lead to the 
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same degree of first-order coherence. Second-order interference effects (and second-order 
cofierence) provide a framework where the statistical properties of the initial state can be 
classified and understood. The final aim of the approach pursued in this paper would be to 
reconstruct, by direct analysis of temperature and polarization correlations, the analog of 
the Mandel parameter and the degree of second-order coherence of the pre-decoupling initial 
conditions. 

To pursue the program briefly outlined in the previous paragraph the first step is to apply 
and translate the theory of (quantum) optical coherence to the case of relic inflaton quanta 
(i.e. relic gravitons and relic curvature phonons). The layout of the paper is therefore the 
following. Section [2] contains a quantum mechanical premise where the second-order correla- 
tions are examined for a single degree of freedom (e.g. a mode of a cavity) but for quantum 
states whose statistical properties are very similar to those arising in the field theoretical dis- 
cussion of the quantized scalar and tensor modes of the geometry. In section E] the quantum 
treatment of the scalar and tensor modes of the geometry is specifically discussed in a unified 
perspective and by emphasizing those aspects which are germane to the present analysis. In 
section H] the degree of first order coherence is computed and analyzed with particular atten- 
tion to wavelengths larger than the Hubble radius. The intensity correlations are studied in 
section El In sections |6] and [7] the degree of second-order coherence is computed in different 
situations and always in the framework of the ACDM scenario. The possibilities of a direct 
estimate of the degree of second-order coherence are also outlined. Section [8] contains the 
concluding remarks and the perspectives of forthcoming analyses. 



By defining the quantum averages with respect to the state |s), the normalized degree of 
second-order coherenc^ of Eq. (11. 7p can be written as 



By coarse graining over technical details which will be the subject of the forthcoming sections, 
it is fair to say that the quantum state of relic phonons (or relic gravitons) belongs to the 
same class of generalized coherent states which arise in the quantum theory of parametric 
amplification of Glauber and Mollow (see, e.g. the Hamiltonian of Eq. (3.3) in the first 
paper of [3Tj). The state \s) introduced in Eq. (12. ip will then be the result of the action 
of a given unitary operator W on a given initial state |in), i.e. |s) = W|in). In this analogy 

For sake of conciseness, the arguments of g^^-* shall be omitted and it will be understood that y*-^-* refers, 
in this section, to a single mode of the field and in the zero time-delay limit. 



2 Single mode of the field 




(2.1) 
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the unitary operator U is to be understood as a time evolution operator. The averages over 
|s) can then be made exphcit. Since = W and [&, = 1, the hnear relation between 
the creation and annihilation operators (a, a)) and (6, 6^) can be parametrized by the two 
complex coefficients c±: 

a = U^bU = c+b + c*^P- (2.2) 

since jc+p — |c_p = 1, c± depend upon three three real numbers (i.e. one amplitude and 
two phases); this occurrence is related to an underlying S'f/(1, 1) dynamical symmetry (see 
e. g. [22] which will be made explicit in section [31 

If the total number of efolds greatly exceeds the maximal number of efolds presently 
accessible by large-scale observation^ (i.e. Ntot ^ A/'max), the state |in) coincides, in practice, 
with the vacuum, i.e. 6|in) = 0. In this case the state \s) preserves the minimum uncertainty 
relations |2Il 1321 ESI El] (see also [3S])- These states are often dubbed squeezed [311E51ES1EZ] 
and lead to a specific degree of second-order coherence which will be extremely relevant 
for the forthcoming considerations. The value of A^max can be computed once the post- 
inflationary thermal history is sufficiently well specified (see, e.g. |5l EH EH]) and will be 
discussed in greater detail later on; for the moment it suffices to posit that A^max — 63. If 
iVtot — A^max — A^min the Statistical properties of |in) can have an impact on the large-scale 
power spectra as argued in [27] and, in this case, |in) does not necessarily coincide with |0). 
Using Eq. (12. 2 p the numerator and denominator of Eq. (12.11) are separately given b}|§: 

{s\a) a) aa\s) = |c_|^(in|666'^6'''|in) + |c4.|^(in|6''" 6''" 66|in) 



I I |2| |2 

+ |c_| |c+| 



{m\bb^ 66^ I in) (in 1 6 6^ b^b\m) 



+ (in 1 6^ 6 6"^ I in) + (in 16^6^1 in) 



s|a^d|s)^ = |c_|^(in|6&^|in) + |c+|^(in|6^ 6|in) 



(2.3) 
(2.4) 



Excluding, for the moment, the logical possibility that the initial fluctuations have nothing 
to do with quantum mechanics (see section [6]), the state |in) can be either pure or mixed. 

Let us now pause for a moment and let us elaborate on the distinction is between pure (or 
mixed) states and correlated states which will be relevant for the forthcoming considerations. 
The purity of a state does not determine, by itself, the degree of second-order coherence. A 
correlated state (i.e. a state exhibiting a degree of second-order coherence potentially larger 
than a coherent state) can be however modeled in terms of a mixed state whose statistical 
weights are appropriately chosen. Consider, for instance, the following parametrization of 



^In the standard terminology it is customary to introduce also A'niin, i-C- the minimal number of efolds 
necessary to fix the problems of the standard big-bang cosmology 
^Units h = c = 1 will be used throughout. 
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the density matrix 

oo 

p= |in)(in| =5]P„|n)(n|, E = 1- (2-5) 

n n=0 

Thanks to the parametrization of Eq. (12.51) . the second-order correlation effects of the state 
I in) will be reflected in the correlation properties of the statistical ensemble defined by the 
weights Pn [IQ]. A non vanishing initial degree of second-order coherence 

(in|6t6|in) ^ ' ^ ^ 

implies, in terms of the parametrization of Eq. (12. 5p . that the statistical weights must satisfy 
the condition 



1 2 



(2.7) 

n=0 '-n=0 -' 

while there are various ways of satisfying the condition (12. 7p it can be shown that if we want 
all the cumulant moments of the distribution Pn to depend only upon the lowest two [H], 
then Pn must obey the following recurrence relation 

(n + l)P„+i = (a + bn)Pn, b ^ 0. (2.8) 

(2) 

If 6 = 0, then 'gU = 1 and P„ is given by the standard form of the Poisson distribution with 
averag^ multiplicity n = a. If, however, 6 7^ the generating function of the distribution is 
simplju 

Ms)-^fP. = j^„ (2.9) 

where clearly = 1 as it must be to be compatible with Eq. (12. 5p . From Eq. (12. 9p 

the various moments of the distribution are obtained by taking the derivatives of Eq. (12. 9p 
at s = 1. It is useful to parametrize the variance in terms of the ratio between a and b 
(i.e. ( = a/b) and in terms of n (the average multiplicity). In terms of these quantities the 
generating function of Eq. (12. 9 p can be written as 

1 1 

•M(^)= r-n ^ = = + 7- (^.lO) 

[n{l — sj + n'^ n Q 

where D"^ and n are defined, respectively, as D"^ = A^"(l)-|-A^'(l) — [A^'(l)]^ and n = A^'(l); 
the prime in the two preceding expressions denotes a derivation with respect to s. While n 



^The probability generating function A4{s) can be directly obtained from Eq. p.Sp even without knowing 
the explicit form of P„ ; it suffices to multiply the right and left hand sides of Eq. (I2.8P by s" and to sum over 
n both sides of the resulting equation. In this way the finite difference equation (|2.8p will be transformed in 
a differential equation in s which can be solved by imposing the boundary condition A4{1) ~ 1. 



11 



simply denotes the number of particles of the initial state, D"^ (and hence Q measures the 
degree of second-order coherence of the initial state. Equation ( I2.10p shows that the higher 
moments of the distribution are all expressible, as anticipated, solely in terms of C and n 
since they are given, by definition, as derivatives of M.{s). The final expression for Eq. (12.11) 
is then given by 



1 rfN +rf{N + if + ANrf{N + 1) 
C [2Nn + n + NY 



N{N + l)[4rl2 + 8n + 1] + N'^{n'^ + 4n + 2) + (A^ + ifri^ 

[2Nn + n + T^f ' ^ ' ' 

where = |c_p and |c+p = 1 + A^. There are various notable limits of Eq. (12. lip : 

• if 7T = Eq. (12. lip reduces to 

^(2) = 3 + 1, (2.12) 
which is the result expected in the case of the squeezed vacuum state; 

• if = and C = 1; '9^'^^ = 2 which is the case of a thermal state (see, e. g. [151 [16]): 

• if A^ = and C oo then 'g^'^^ = 1, which is the case of a coherent state (see, e.g. 

mm); 

• if iV > 1 and n > 1 then ^(^^ 3(1 + l/C)/2; thus, if C -> 1 we shall have that also 
^(2) ^ 3 as long as A^ ^ 1 and n ^ 1: this result agrees with the well known result 
concerning the (single mode) squeezed thermal states ^2] . 



The most general situation corresponds to the case where not only n ^ and A^ 7^ but 
also when C 7^ 0. It is appropriate to conclude this discussion with a simple remark on 
the normal ordering in the definition of the degree of second-order coherence. Suppose that 
we define the degree of second-order coherence without resorting to normal ordering; for 
instance we can antinormal order (i.e. (s|dda^ (^H^)) o^' even adopt a mixed kind of ordering 
(e.g. (s|d d^ d^ d|s)). Denoting with g^^^ the degree of second-order coherence with generic 
ordering, it can be shown that 

^(2)=^(2) + o(i/(Ar)), (2.13) 

where (A^) schematically denotes the mean number of quanta which depends, ultimately, 
upon the initial state. 
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3 Relic phonons and relic gravitons 



The simplest setup compatible with the ACDM paradigm stipulates that the background ge- 
ometry is conformally flat with metric tensor g^,^ = a'^(T)r]^i,, where r denotes the conformal 
time coordinat^. The scalar and the tensor fluctuations of the geometry can be described 

asm mm 



^{s)goo = 2a^0, S(^s)gij = '^a^'ip ^i. 



S{t)9ij = -a hij, d,h] = h] = 0, (3.1) 

where S(s) and 6(t) denote, respectively, the scalar and the the tensor fluctuations of the 
corresponding quantity; the gauge freedom has been completely flxed in Eq. (13. ip by selecting 
the conformally Newtonian gauge. If the inflationary stage of expansion is driven by a single 
background scalar fleld ip, deflning with 5(s)V5 the scalar fluctuation of the inflaton and with 
ip the actions describing the evolution of the scalar and tensor modes of the geometry can 
be written, respectively, as (see, for instance, flO\ [TT]) 

5(3) = \j d^x ^^r^d^n dpU, (3.2) 
'^'(t) = ^ j d'^x ^f^^g^'^dahijdphij, (3.3) 

where TZ represents the curvature perturbation on comoving orthogonal hypersurfaces whose 
explicit expression, in the case at hand, is 

Orif Or In a 

and T-i = 9^ In a. Using that g^,^ = a?{T)ri^y, Eqs. ( 13. 2p and (13. 3 p can be written as 

5m = \j d^x t^dJldpUz^ (3.5) 



S(t) 



^Jd'^xr]"^do,hd^ha\ (3.6) 



Equation (13. 6 p holds for each of the two tensor polarizations (see, e. g. [101 SS]) having 
deflned = v^£p/i with ^ = (B, ® and ip = l/VSirG (see also Eq. (16. 5 p where the reduced 
Planck mass Mp = ip^ enters the deflnition of the slow-roll parameters). Indeed hij{x,T) 
can be decomposed as 

hij{x, t) = J2 Qifhxix, r), (3.7) 



^The conventions adopted in the redshifts are such that the present value of the scale factor oo is nor- 
malized to 1; note that in section [T]t denoted, consistently, the time coordinate in Minkowski space 
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where g® = (ajSj — bibj) and qfj = {hibj + aj&j); defining with k the direction of propagation 
of the wave, a, h and k form a triplet of mutually orthogonal unit vectors. Probably the first 
paper mentioning quantum mechanics as a possible source of large-scale inhomogeneities, 
though not in the framework of any infiationary hypothesis, is the one of Sakharov [47] . 
The emphasis on the action for the normal mode of the scalar fiuctuations (i.e. scalar 
phonons) appeared in a paper by Lukash [48j in the context of fiuid models. Later on different 
authors applied it to scalar field matter with particular attention to the quantization of the 
fiuctuations [19], |50] (see also [51]). The form of the actions is the one derived in [10]. The 
scalar and tensor fiuctuations of the geometry can be canonically quantized; after introducing 
the appropriate normal modes 



/i(x, r) = a(r)/i(x, r), z/(x, r) = z{t)71{x, r). 



the tensor and scalar Lagrangian densities become, respectively, 

£(t)(f, r) = i {drfif + {dr In a)V' - 2(5, In a) - {d.fif 
£(^)(f,r) = i {druf + {dr\nzyn^ -2{drlnz)udriy-{diuf 

whose associated canonical momenta 



7r(t) = f^r/U — (9, In a)/i, 7r(s) = dri^ — {drlnz)^, 



can be used to derive the canonical Hamiltonians 

1 



H(t) = 2 / ^^^ + 2(5r lna)7r(t)/i + {difif 
H{s) = ^ / ^^^ ^?s) + 2(<9, lnz)7r(s)Z/ + {diuf 



(3.8) 

(3.9) 
(3.10) 

(3.11) 

(3.12) 
(3.13) 



Since Eqs. f l3.12p and f l3.13p have the same canonical structure, the two problems can be 
treated simultaneously by resorting to the Hamiltonian: 



H(t) 



(TX 



(3.14) 



where $ and ii are the two canonically conjugate field operators. In Eq. fl3.14p A = 
i{dr\n.a) /2 (in the case of the tensor modes) and A = ^((9,- In z)/2 (in the case of the scalar 
modes). Similarly tt will coincide either with 7r(t) (in the case of the tensor modes) or with 
7r(s) (in the case of the scalar modes). The Fourier representation of the field operators can 
be written as 



$(f, t) 



1 



r e 



-ip-x 



1 



r e 



-ip-x 



(3.15) 
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where V represents a fiducial (normalization) volume. In the continuum limit we will have 
V J d^k/{2TTY and the canonical commutation relations impose, in Fourier space, 
[•f*^, vrt] = iS^^\k — p) where $t = <l'_^ and ttI = 7r_^ because of the hermiticity of the 
corresponding field operators in real space. Introducing creation and annihilation operators 
obeying [a^, at] = S^^^k—p), the field operators and the canonical momenta can be expressed 



as 



(3.16) 



Inserting Eq. fl3.16p into Eq. fl3.14p the resulting Hamiltonian in the continuum limit is 



A*(r)/C_(p-) + A(r)/C+(p-) 



(3.17) 



^o(p) = 2 



(3.18) 



Hir) = 2 j d'p^ p ICoip) + 

where the operators }C±{p) and /Co(p) 

JC+{p) = ata^., IC-{p) = a^a_p-, 

satisfy the commutation relations of the SU{1, 1) Lie algebra, i.e. 

[JC4p), JC4q}] = 2 /Co(p) 6^'\p- q), [/Co(p1, /C±(g-)] = ± IC^^p) 5^'\p - q). (3.19) 

The group SU{1, 1) is not a symmetry group of the Hamiltonian of the problem but the 
SU{1, 1) algebra can be viewed, in the terminology of [52] (see also [53]), as the spectrum 
generating algebra insofar as the total (generalized) charge does commute with all the gener- 
ators of the group (as well as with the total Hamiltonian) while the total number of particles 
does commute with the charge but not with the full Hamiltonian. Owing to the group struc- 
ture f l3.19p and to the specific form of the Hamiltonian of Eq. ( I3.17p . the multiparticle final 
state can be obtained by applying to the initial state \^i{p)) the product of two unitary 
operators S((y9p) and S(crp): 



\mf{p)) = S(^,) S(a,)|vl/,(p-)), |vl/^) = n |^/(P)), 



(3.20) 



where the unitary operators are defined aqj 

n{(pp) = exp[-2i(ppKo{p)], J:{ap) = exp[a*K_{p) - (TpK+ip)], (3.21) 

with cTp = rpC^^p and ap = {2ipp — 7p); the time evolution of the variables rp(r), ^Ppij) and 
Oipij) is given by 

— — = ZiA cosa„, 
cLt 



^Note that tpp should not be confused with tp (denotmg the inflaton field). This confusion cannot actually 
arise since it is clear that (pp is a momentum-dependent phase. 
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difp ^ . , , 

— = p — 2iX tann rp sm «„, 
dr 

-^ = 2p-4zA— (3.22) 
dr tann 2rp 

Equations fl3.22p are symmetric for — )■ — and A — t- — A. The corresponding Hamiltonian 
of Eqs. f l3.14p is actually symmetric for 

and analogously in the tensor case for a — )■ 1/a. The transformation of Eq. fl3.23p is related 
to electric- magnetic duality [HU [551 iSE] in conformally flat background geometries when the a 
and z are replaced by the (dynamical) gauge coupling. The relation between the Schrodinger 
and the Heisenberg descriptions is easily worked out by appreciating that 

H^) = fk{T)d^ki^,) + fldUir,). (3.24) 
7rfc(r) = c/fc(r)a^(ro) + ^fca^_g(ro). (3.25) 

In section |2] the operators at the initial time tq have been denoted with 6; thus, accounting 
for the momentum dependence, we will denote 6g = a^(To), 6^- = a^^(ro). The evolution 
of the mode functions /^(t) and gk{T) can be obtained from the evolution equations in the 
Heisenberg description: 

idr^ = [$, H], tdrir = [tt, H], (3.26) 

and they are given by 

drfk = gk + dr In afk, d^Qk = -k'^fk - {dr In a)gk (3.27) 
for the tensor case and by 

drfk = gk + dr In zfk, drh = -k'^fk - {dr In z)gk (3.28) 

for the scalar case. Both sets of mode functions are subjected to the following Wronskian 
normalization for any r 

fpir)~g;ir)-f;iT)~g,iT) = t, (3.29) 
fpir)g;iT) ~ f;{T)g,iT) = z. (3.30) 



4 First-order coherence 

The degree of first-order coherence is measured by the following normalized correlation func- 
tion: 

g^^\.^ylr)^^mMl]^. (4.1) 
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In the usual quantum optical setup provided by standard interferometers (such as Michelson, 
Mach-Zehnder or Sagnac) the ^-operators are replaced either by electric field operators or 
even by classical fields. In Young interferometry the electric fields are split, time delayed and 
then recombined on a screen. It is easy to show that the analog of Eq. (14. ip is nothing but 
the visibility, i.e. the normalized difference between the maximal and the minimal intensity 
of the light detected on the Young screen. 

In quantum theory, the degree of first-order coherence is fully determined by the averaged 
multiplicity of the initial state. Let us suppose that the initial state is characterized by an 
averaged multiplicity Wq per each field mode. Following the notations employed in section [2] 
but accounting for the momentum dependence we will have 



{m\bibj^\m) 



n, 



(4.2) 



Given the average occupation number of the initial state, Eq. ( 14. ip can be computed in 
explicit terms. Recalling Eqs. ( I3.15P and ( I3.16P the two point function is simply given by 



(l>(f,r)<l(f,r)) 
(: <l(x,r)l'(f,r) :) 



47r2 



q dq 



cosh 2rq — cos aq sinh 2rq 



(2rT, + l)jo(gr), (4.3) 



27r2 



q dq 



cosh 2rq — cos aq sinh 2rq 



(4.4) 



where the normal ordered case has been included for comparison and where jo{qr) is the 
zeroth-order spherical Bessel function [571 [58]. To deduce Eqs. (14. 3 p and (14.40 it is useful to 
recall that 

(4.5) 



aq = e 



cosh Tq bq' — e*"^' sinh rqb^_^ 



where, as in Eqs. (I3.20p -( 1X^ . 7^ = {2(pq — aq). Since sinh Vq 



Nq is the averaged 



multiplicity of the squeezed vacuum state, the expression appearing inside the square brackets 
in Eqs. (14. 3 p and (14. 4p can also be written as 



cosh 2rq — cos aq sinh 2rq = 2Nq + 1 — 2y Nq{Nq + 1) cos aq. 



(4.6) 



The solution of either Eqs. ( I3.22p or Eqs. (I3.27p -( l3.28p leads to an even more explicit form 
of the degree of first-order coherence which will be indirectly mentioned in section [7[ In 
spite of the statistical properties of the initial state and in spite of the operator ordering 
Eqs. ([13]) and ([13]) imply that 



lim g^^\r, r) 

qi — s-0 



F - y\ 



(4.7) 



Concerning the limit of Eq. (14. 7p few comments are in order. Mathematically the correct 
limit to be implemented is r — )■ since the correlation function is the result of an integral 
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over the comoving three- momentum; at the same time the physical limit, as indicated in Eq. 
(14. 7p is gr ^ 1 (and also, as we shall see in the case of second-order correlations, |A;r| <^ 1). 
From the explicit expression of g^^\r,T) it is also clear that the integrals over q might not 
always be convergent. Still the result of the limit holds because, for r — )■ the divergent 
contributions in the numerator and in the denominator exactly cancel. 

The result of Eq. (14. 7p can be dubbed by saying that the relic gravitons and the relic 
phonons are always first-order coherent irrespective of the statistical properties of the 
initial state which could be rather different such as a mixed or a pure state. In both cases 
the density operator of the initial state can be defined in the most appropriate basis, for 
instance a Fock basis or a coherent state basis. For a thermal (or chaotic) ensemble the 
density matrix can be written, in the Fock basis, as 

— "fc 

P = EA"}IM)(WI' P{n} = I[7T-^y^, (4.8) 

where, in analogy with the notations employed in section |2l = Tr[patdg] is the average 
occupation number of each Fourier mode and, following the standard notation, \{n}) = 
l^fei) ll'^A;2'^ W^k^)--- where the ellipses stand for all the occupied modes of the field. The 
density matrix always describes a mixed state but the should not be necessarily identified 
with the Bose-Einstein occupation number. In the case of a multimode coherent state the 
density matrix can instead be written as 

k 

m}) = e-l^^^'^/^E^I^A?)- (4-9) 

The initial state of Eq. (14. 9 p is pure and its averaged multiplicity per Fourier mode is given 
by rlfc = Hence Eq. (14. 7p holds both for the state of Eq. (14. 8 p and for the state of 

Eq. (14. 9 p as well as for all the possible initial states (either pure or mixed). This simply 
means that the degree of first-order coherence is only sensitive to the particle content of the 
initial state but not to its statistical properties. Initial states exhibiting a high degree of 
correlation cannot be distinguished just by looking at the analog of Young interferometry. 
A given multiparticle state can always be projected on the coherent state basis [15]. 

^''The terminology "first-order coherence" (and later on of "second-order coherence") is the one borrowed 
from the Glauber theory of optical coherence as formulated in [25j [26] . In section [T] the basics of Glauber 
approach to optical fields have been reviewed and will now be applied in the present and in the following 
sections to the case of relic phonons and relic gravitons. 
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The idea is to write the multiparticle density matrix of a mixed state in the basis of a 
multimode coherent state using the overcompleteness of |{7}) 

P = I P{{i])\{i]){b]W{i] (4.10) 



dpii} = n 



(4.11) 



where P{{'~f}) is the phase-space functional. The representation of Eq. fl4.10p cannot be 
more singular than a Dirac delta function and it should also be positive semi-definite: these 
two properties are not satisfied by any quantum state. For instance, in the case of the density 
matrices of Eqs. (14. 8 P and (14. 9 P the corresponding phase-space functionals are: 

nil}) = e^p[-i7dv%], (4.12) 

^({7}) = U^^'Kl—P^)- (4.13) 

k 

There are indeed states leading to a phase-space functional which is more singular than a 
delta function (the delta function case corresponding to a coherent state) such as the squeezed 
states and various of their generalizations [59| [60] . As it will be clear in what follows, 
the method of the phase-space functional will just be mentioned as a cross-check in the 
calculation of some expectation values involving the initial state of the relic phonons and of 
the relic gravitons. In the cases where this technique will be employed, the P-representation 
will always be well defined and regular. When normal ordering is imposed, in the coherent 
state basis defined by Eqs. (I4.10p - (l4.1ip and by Eqs. (I4.12p - (l4.13p . the quantum averages are 
replaced by averages over complex numbers weighted by the phase space functional related 
to the so-called Glauber- Sudarshan P-representation [151 [161 [T7] . The latter properties goes 
under the name of optical equivalence theorem and greatly simplifies the calculations of field 
correlators provided the density matrix of the quantum states involved in the average can be 
represented in the coherent state basis with P-distribution not more singular than a Dirac 
delta function [151 [HI [HI- This observation can be used to check the results of various 
correlators in the limit of large occupation numbers since, in this limit, the ordering of the 
higher-order correlators is immaterial as discussed at the end of section [2 

The main result of this section can be summarized by saying that the field describing 
the relic phonons and the relic gravitons is always first-order coherent when the relevant 
wavelengths are larger than the Hubble radius at each corresponding epoch. Furthermore, 
given the explicit form of Eq. (14. 3 p and (14. 4p . we can also conclude that < g^^\r, r) < 1, 
i.e. the field is first-order coherent in the large-scale limit and partially coherent for smaller 
scales. 



19 



5 Intensity correlations 



The degree of second-order coherence is defined as in Eq. (11.61) . i.e. 

(2)/-- X (X(x, r) X(y,r)) 
g^'[x,y;T) = — , (5.1) 

where X(x, r) = $^(x, r) and X{y, r) = r). As anticipated in the introductory section 
Eq. (15. ip differs from the quantum degree of second-order coherence appearing in Refs. 
[13 [151 El- In quantum optics the intensities are measured by phototubes and the correla- 
tion is proportional to the transition rate for a joint absorption of photons at the two points. 
The treatment of the photoelectric effect shows that the transition amplitude is proportional 
to the matrix element of E^^HiJ, r?) E^^Hx, ri); accordingly the degree of second-order co- 
herence is defined, in quantum optics, ag^^l[TS| [TB| [T7] 

' " (:^(-)(f,ri)i(+)(f,r0:)(:^(-)(y,r2)i(+)(y,r2):)' 

where E^~\x, r) and E^~^\x, r) denote, respectively, the negative and the positive frequency 
parts of the electric field operator for a single polarization. By rewriting Eq. (15. 2 p in the 
notation of the present paper the quantum degree of second-order coherence becomes, in the 
case Ti = T2 = T, 

. (: l>(")(f,r) l>(-)(f,r)l>(+)(^,r)l>(+)(x,r) :) 

'[x,y\T] = , 5.0 

(: $(-)(f,r)$(+)(f,r) :)(: $(-)(y,r)$(+)(y,r) :) 

where, as above, <l>''^'*(x, r) and (x, r) denote the negative and the positive frequency 
parts of $. Equations (15. ip and (15. 3p are technically different but physically equivalent. 
The numerical value of the normalized degree of second-order coherence (in the zero time- 
delay limit) can differ between Eqs. (15. ip and (15. 3p for a given quantum state. However, as 
discussed in Eq. (I2.13P these difference vanish either when the number of particles of the 
initial state is large or when the number of produced particles is large. Barring for specific 
numerical differences which are relevant in the limit of small occupation numbers, Eq. (15. Sp 
shall be primarily considered; if appropriate, the relations of the obtained results with the 
normal ordered definition shall be swiftly mentioned. The degree of second-order coherence 
given in Eq. (15. ip can be estimated as 

(X(x, r)X(y, r)) = E E E E J'i^^ V. q' , p'; ^, r) (5.4) 



^-'^As in section [21 g^^^ denote the normal ordered degree of second-order coherence while in Eq. ()5.ip the 
bar has been omitted. 
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where the expression p, g', p'; x, y, r) is given by 

{ala^a^. at, )e^(9--P)---i('r'-p')-y + (aj-ap^at, at, )e^(9-P)--+i('r'-P')-j/ 1. (5.5) 



Inserting Eq. fl4.5p into Eq. (15.51) various averages will appear such as 

(hVhlk\^,), (5.6) 

where the ellipses stand for the four remaining permutations. Each o the averages Eq. (15. 6p 
are evaluated using the density matrix of the initial state, i.e. for instance 



blb^' bp') = Tr 



phVhlh^.hj;' 



(5.7) 



and similarly for all the other expectation values of the fields arising in Eq. (15. 5p upon 
insertion of Eq. (14. 5p . The density matrix appearing in Eq. (15. 7p can be the density matrix 
either of a pure state or of a mixed state and can be written, in general terms, as 

P = E An}IW)(WI, EAn} = l- (5.8) 

{n} {n} 

As already mentioned in section [2] we are interested in the possibility that the initial state 
has a specified degree of second-order coherence. The simplest non-trivial situation is the 
one discussed in Eq. (12. 8 p and hereby generalized to the situation where a — a^ and 6 — )■ fefc 
do depend upon the momentum (as appropriate in the case of many bosonic degrees of 
freedom). Thus the field theoretical generalization of the probability distribution implicitly 
mentioned in section [21 i.e. 

^"^"Vr(c,)r(n, + i)U, + C,J U, + C,J • ^'-'^ 

The evaluation of Eq. (I5.7p proceeds therefore by noticing that 



{h\ h] bk hi) = {hi bj k bi)6ij 5jk5ik + 

{h\ h] k k) ^rk^ji[l - 5.,] + {h\ h] b, h) S^eSjk[l - S,,], (5.10) 



where hi and 6] denote the annihilation and creation operators related to two generic mo- 
menta, i.e. for instance 6^- and fet; furthermore, following the same shorthand notation, Sij 
denotes the delta functions over the three-momenta (i.e. 5g,p). In Eq. (I5.10p two different 
classes of terms appear: in the first class of terms (i.e. lj)\h\hihi)) there are four different 
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operators all acting on the same momentum; in the second class of terms (i.e. {b] hi bj) and 
(bj bj bj bi)) the momenta are paired two by two (as the corresponding deltas indicate). Since 
the averages are to be computed using the initial density matrix characterized, in general, 
by the statistical weights discussed above we have 

{b\ b] bk bi) = - l)Pi{'mi) 

rrii 

(H 6] hi hj) = Pi j i^i ' "^i ) ' 

{h\h]hjhi) = Y mi rrij Pij{mi,mj). (5.11) 

In Eq. (15.111) the following shorthand notations have been used 

Pi{mi) = PtXrrikJ, PijirUi, mj) = PuX^k,) Pk,{mk^)- (5.12) 

The second relation of Eq. (I5.12p is indeed trivial in the light of the very definition of P{m}'- 

P{m} = Pkiimk^) Pk2{mk^) Pksimks) ■■■■] (5.13) 

where the ellipses stand for the product over the various momenta; at the same time the 
explicit appearance of Pij is useful to trace the origin of the various terms. Inserting Eq. 
(15. lip into Eq. (I5.10p the correlator becomes: 



{hi h] hk hi) = Y mi{mi - l)Pi{mi) - 2mj ^ rrijPij^mi, nij] 

rrij 

SikSji + Sie6jk ■ 

Finally, using Eqs. (I5.12p and (I5.13p . Eq. (I5.14p can be written as 



i Sjk 5 Ik 



+ Y rrii rrij Pi j{mi, rrij] 

nii, nij 



(5.14) 



{hi h\ hk hi) = mii{ini - l)Pi(mj) - 2iniPi{mi) Y ^jPj 



Sij SjkSik 



+ Y^iPii^i)^^jPi 



^ik^jl + SiiSjk 



(5.15) 



For an explicit evaluation of the sums of Eq. (15.150 it is useful to employ the probability 
generating function and the cumulant generating function 



M = n-^A;(^fc'%'Cfc), c = H'^fcl^fc'^'Cfc), 

k k 



(5.16) 



whose specific form, from Eq. (15.90 . becomes: 

Mkisk,n^Xk) 



^k 



lCk + i^-^k)nkr^'- 



^ki^k^n^Xk) = -Cfc In 
22 



nr 



1+{1-Sk)-^ 

'^k 



(5.17) 



and coincides with the expression aheady derived in Eq. fl2.10p in the case of a single degree 
of freedom. The various sums of Eq. flS.lSp can be exphcitly evaluated as combinations of 
the derivatives of the probability generating function. The final result can be expressed as 



(6| bk h) = 72- - 1^6ij 6jk Ski + rii rij 



'jk 



(5.18) 



By restoring, in Eq. (15.181) . the standard notations for the comoving three- momenta of the 
field we will have that: 



,q' 5q',p' +nan 



q 



q,p'^p,q' 



(5.19) 



The same procedure leading to Eqs. (15.181) and (15.191) can be used to compute, with the 
due differences, all the other expectation values arising in the evaluation of the intensity 
correlation. In particular 



{kiybkbi) 

{bi b, bi k) 
{bihybkbi) 
{kb]bih) 



= (^^ ~ 1 j ^ij ^kl + {rii + 1) {Uj + 1) 

+ nk{ni + l)5it5jk, 



n^[-^ - 1 )6ij 6jk6ki + iiiUkSijSke + iii^Uk + l)5ii5jk, 



- l]6ij 5jk6ki + ni{ni + 1) 



^ij^ke + ^ik^e 



- l]6ij 5jkSk£ + niiui + 1) SijSki + SikSij 



(5.20) 

(5.21) 
(5.22) 

(5.23) 

(5.24) 



The physical role of the quantum correlations can be neatly understood by taking the Bose- 
Einstein limit in the correlators of Eqs. (I5.19P and (I5.20p - (l5.24p . In the latter limit Q ^ 1 
for every mode of the field; Eq. (15. 9p turns into the standard thermal ensemble and the 
contributions of the terms of the type (foj b\ b^ bi) exactly cancels; the standard rules of 
evaluating correlators in thermal field theory is quickly recovered [HIl E2|- It is useful, 
at this level, to take the limit Qk ^ oo uniformly for all modes of the field. From Eq. (I5.17P 
it can be immediately appreciated that, in the limit Cfe ~^ oo, the probability generating 
function becomes Aik exp [{sk — l)^fc] which is exactly the generating function of the 
Poisson distribution. The Poisson distribution for each mode of the field is customarily 
associated with a multimode coherent state but this is not exactly our case: the situation 
described by Eq. (15. 8p and (15. 9 p in the limit Cfc — ^ C)0 is the one of a mixed state with 
Poissonian distribution and this cannot be identified with a multimode coherent state since, 
in the latter case, the off-diagonal elements of the density matrix do not vanish (as it is 
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the case for Eq. fl5.9p ). Finally, if the limit (^^ — )■ is taken uniformly for all modes of the 
field, the probability distribution for each /c-mode becomes logarithmic as it can be shown 
by using directly the recurrence relation characterizing the distribution of Eq. (15. 9p . 

We are then in condition of computing explicitly the intensity correlations in terms of an 
initial state characterized by the presence of quantum correlations reducing, in appropriate 
limits, to various statistical mixtures. The result can be written as 



(5.25) 



(X(f,r)X(f +f,r)) = J dink gi^\k,T) [2 + cos krjo{kr)] 

+ Jdlnkgi^\k,T)[l + 2Ukr% 



where jo{kr) is the spherical Bessel function of zeroth order; the two functions Q^^^k, r) and 
^|^)(A;,r) denote, respectively, the bulk (volume) and boundary (surface) contributions 

n 2 



with r 



\x 



Gf\k,r) 
- y\ and 



k nlir) 



d^q 



1 

k^ 



647r5 q\k-q\ 



2Nk{r) + 1 - 2JNk{Nk + 1) cos ak 



^(g,r)^(|A;-g1,r) 



(5.26) 
(5.27) 



^(g, r) = [2 rl,(r) + 1] |2iV,(r) + 1 - 2^Ng{T) [Ngir) + 1] cosa,(r)|. (5.28) 

We shall also assume, from now on, that = Uk and that Q = (k with \k\ = k. The 
denominator appearing in the degree of second-order coherence is 



(X(x,r))(X(y,r)) 



d\nkg^^\k, 



47r2 



(2nfc + 1) 



2Nk + l- 2dNk{Nk + 1) cosafc 



(5.29) 
(5.30) 



where the superscript reminds that we are dealing here with the square of the first-order 
correlation discussed in section [31 The normally ordered definition of the degree of second- 
order coherence introduced in Eq. (15. 3p leads to a result which is somehow similar to the 
one obtained in the case of Eq. (15. ip . For future comparison the normal ordered intensity 
correlation can be written, in explicit terms, as 



(: l>(-)(x,r)<|)(-)(y,r)<|)W(y,r)<|)(+)(f,r) :) 



+ 



rflnA;^;^ (A;,r) 

d\nkgf\k,T)[l+3o{kr)l (5.31) 



where now 



knl n 



47r2FVa 



2Nk + l- 2JNk{Nk + 1) costtfc 



(5.32) 
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647r5 q\k-q\ 



2Ng+l- 2jNg{Ng + l) 003 0, 



(5.33) 
(5.34) 



Of course in case the normal ordered expression is used, also the denominator of Eq. (I5.3p 
must be ordered and computed accordingly. The result is 

2 



d\nkg^^\k,T 



2n' 



nk 



2Nk + l-2JNk{Nk + l) cos«fc 



(5.35) 
(5.36) 



There are two main differences between the normal ordered correlators and the non-normal 
ordered ones: a numerical factor in Eq. fl5.3ip and the ubiquitous presence of Hk instead of 
{2nk + 1). The expressions of Eqs. fl5.32p and f l5.33p vanish in the limit — )■ while in the 
case of Eqs. (15.261) and (15.271) the same limit does not vanish. Furthermore, the degree of 
second-order coherence inherits a volume dependence which is directly linked to the existence 
of an initial state with a non vanishing degree of second-order correlation. To conclude the 
present section it is appropriate to mention that the averages over the initial state can also 
be conducted by making appropriate use of the phase-space functional previously discussed 
(see Eqs. fl4.10l) - fl4.11l) ). In particular, the multiparticle states defined by the density matrix 
of Eqs. (15. 8p and (15. 9p leads to a P-representation given by 



^({7}) = n 



|^^|-2(l-Cfc) g-|7fcPCfe/nfc_ 



(5.37) 



The P-representation of Eq. (I5.37P has interesting limits. In particular for (^^ — )■ 1 the P- 
representation of Eq. (14.120 is reproduced. The representation (15.370 has been first discussed 
in [63] in the purely quantum mechanical case. 



6 Degree of second-order coherence 
6.1 Basic considerations 

The degree of second-order coherence discussed in section O will now be evaluated explicitly in 
various potentially interesting situations. The specific values of the cosmological parameters 
determined using the WMAP 7yr data alone in the light of the vanilla ACDM scenario [U 
are0: 

(fib, f^c^de, /io, ^rc) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088), (6.1) 

^^Note that Ere denotes the optical depth at reionization and has nothmg to do with the slow-roll parameter 
e which will be introduced in a moment. 
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where Qx denotes the present critical fraction of the corresponding species (i.e., respectively, 
baryons, CDM particles, dark energy); Ug denotes the scalar spectral index (see also Eqs. 
(16. 4p . (16. 5 p and (16. 6p hereunder). If a tensor component is allowed in the analysis of the 
WMAP 7yr data alone the relevant cosmological parameters are determined to be [HE] 

(fib, ficfide, ho, ns. Ere) = (0.0430, 0.200, 0.757, 0.735, 0.982, 0.091). (6.2) 

In the case of Eq. (16. ip the amplitude of the scalar modes is An = (2.43 ± 0.11) x 10~^ 
while in the case of Eq. (16. 2 p the corresponding values of A-jz and of are given by 

= (2.28 ±0.15) X 10-^ rt<0.36 (6.3) 

to 95 % confidence level. The experimental parametrization of the scalar and tensor power 
spectra is [H [2] 

Vn{k) = An[^J , n(fc)=A(^) , (6.4) 

where rit, are, respectively, the tensor spectral index, and kp = 0.002 Mpc~^ is the pivot 
scale; = At /An denotes the ratio between the tensor and the scalar power spectrum 
at kp. The qualitative features of the effects discussed here do not change if, for instance, 
one would endorse the parameters drawn from the minimal tensor extension of the ACDM 
paradigm and compared not to the WMAP 7yr data release but rather with the WMAP 3yr 
data release [TH [72], implying, for instance. An = 2. ills ^ 10~^, Ug = 0.984 and rt < 0.65 
(95 % confidence level). 

The numeric values reported in Eqs. (16. ip . (16.20 and (16.30 determine bounds on the slow 
roll parameters appearing directly in the evaluation of the degree of second-order coherence. 
In particular, within the present notations, the slow-roll parameters are defined as where, as 
usual, 



H"^ 2 \ V 

7] = ^ = €-7], rj = Mp-j^, (6.5) 

and Mp = Mp / y/Sn is the reduced Planck mass, H is the Hubble parameter, V is the 
inflaton potential and the overdot denotes a derivation with respect to the cosmic timt^ 
coordinate t. To lowest order in the slow-roll expansion we also have that 

rit = — 2e, Us = 1 — 6e + 2r/, rt = 16e = —Srit. (6.6) 



^■^Recall that, as usual, the relation between cosmic and conformal time parametrization is given by 
a(T) dr — dt. 
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The mean number of particles per Fourier mode can be obtained from Eqs. fl3.22p . While 
Hk and (k are a property of the initial state, the quantities Nk and at can be computed 
explicitly, for instance, in the case of a single-field inflationary model. To get the explicit 
solutions either in terms of Eqs. (13.220 or in terms of Eqs. (I3.27p - (]3.28p background fields 
need to be expressed in terms of the conformal time coordinate and in terms of the slow-roll 
parameters of Eq. (16. 5p : 



d-r In a 

d^^\na+{drlnaf 
d^lnz + idrln zf 



aH 



1 



r(l-e) 
l + e + 7] 



r(l-e)' 
2-e 



2 + 2e + 3r] + er] + r]^ 



(1 



(6.7) 
(6.8) 
(6.9) 
(6.10) 



Recalling that A^^ = sinh^ and denoting with A^^^^ and A^^*^^ the average number of phonons 
and gravitons the final result turns out to be: 



(1) 



\Hj^\-kr)f + \H;^\i-kr)\ 



(6.11) 
(6.12) 



where if^^^ (z) is the Hankel function of the first kind of index 7 and argument z [STJ [58] . 
In Eqs. (16. lip and (I6.12p the indexes of the Hankel functions depend upon the slow-roll 
parameters as 

P-'-^^. (6.13) 



2 1 



2 1 



Using Eq. ( 16. 6p a and (3 can be written directly in terms of the scalar and tensor spectral 
indices 

6 + nt „ 8 + 3nt — 2ns 



a 



(3 



(6.14) 



2(2 + nO' " 2(2 + nO 
In the range of parameters mentioned in Eqs. (16. ip . (16.20 and (16.30 the indexes a > 1 and 
/3 > 1 are always positive definite and, moreover. 



(3 - 2a) 
(3 - 4a) 



2e 

"1 - e 
3 + e 
1 -e 



<0, 
<0, 



(3-2/3) 
(3-4/3) 



4e + 21] 

3 + 9e + 4r7 
fl-e) ■ 



(6.15) 



Equations (16.110 and (I6.12p can be expanded in the limits |A;r| ^ 1 and |A;r| ^ 1 holding, 
respectively, when the corresponding wavelengths are either shorter or larger than the Hubble 
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radius. In the limit |A;r| <^ 1, Eqs. fl6.11l) and fl6.12l) become: 



2n 
2tt 



1-2/3 
^^1-2/3 



1 + 



4(/3-l)2 
|A;rP 



Ma - 1] 



(6.16) 
(6.17) 



In the limit of short wavelengths the particles are all inside the Hubble radius and, conse- 
quently, 



(p) 



(6.18) 



where the factor 1/2 confirms, a posteriori, the correctness of all the normalizations and 
implies that, in the vacuum, there is half a quantum per Fourier mode. The degree of 
second-order coherence is also determined by a phase (see, e.g. Eq. ( I5.28P ) whose explicit 
form for phonons and gravitons can be written as 



cos a 



(p) 



vr, , \Hj^\{-kr)f - \Hj^\-kT)\' 

-{-kT}^^ 



COS a 



(g) 



TT 



\H. 



ivf)(ivf) + i) 



(6.19) 
(6.20) 



where, consistently with the notation of Eqs. (16. lip and (I6.12p . the superscripts refer, 
respectively, to the case of the gravitons and of the phonons. In the limits \kT\ <^ 1 and 
|A;r| > 1, Eqs. fl6J9D and flOOD lead, respectively, to 



lim cos a 

\kT\->-0 



(s,t) 



^ -1, 



lim cosa[.'^'*'' — )■ 

|fcT|— >00 



(6.2i; 



where the superscripts simply mean that the mentioned limits hold separately for the scalar 
and for the tensor modes. 



6.2 Initial vacuum state 

The relevant physical limit of the degree of second-order coherence can be expressed as 
follows: 

lim q^'^Ux,x + r:T). (6.22) 

The limits appearing in Eq. fl6.22p imply that the degree of second-order coherence is 
evaluated at coincidental spatial points and when the comoving momenta are much shorter 
either than dr In a (in the case of the gravitons) or than dr In z (in the case of the phonons). 
As in the case of the limit discussed in Eq. (14.70 the mathematical definition of the limit 
only involves r since the degree of second-order coherence is always defined as the ratio of 
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integrals over the comoving three-momentum. At the same time the limit r — )■ is physically 
realized for typical wavenumbers \kr\ — )■ and |A;r| — )■ and this defines, according to Eq. 
f l6.22p the normalized degree of second-order coherence in the asymptotic limit of typical 
scales (or wavelengths) larger than the Hubble radius. This kind of procedure has some 
analogy with quantum optics where often the results of the momentum integrations in the 
numerator and in the denominator of the degree of second-order coherence simplify in the 
limit r — )■ and lead to a degree of second-order coherence which only depends on time (see, 

e.g., mm- 

It is useful to remark that, a posteriori, the order of the limits appearing Eq. (16.221) is 
not essential: the same results can be obtained by changing the order in which the limits 
are taken. From the explicit expressions of Eq. (15.251) the degree of second-order coherence 
can be written as 

recalling that jo{z) 1 for z ^ 0. Using all the results derived in the present section (and, 
in particular, Eqs. (16.161) . (I6.17P and (16.211) ). also the limit |A;r| — )■ can be taken with the 
result thatl3 

lim g^'^\r,T) = A(a) + 'Bs(a), (6.24) 

I A^T I " — ^0 J I /t^ 7" I — ^0 

A (a) = ^/^^MMil^M^i^ .6 25) 

' V |/dlnA;(2rrfe + l)A;3-2«|' ' ^ ' ^ 

^ ^ 3 jd\nkejd\n qq^-'- (2n, + l)y{\k - q]) 
'^"^ 2 |/c/lnA;(2nfc + l)A;3-2"|2 > I • J 

where _ 

y(^\k-q\)= dx^^ , \k-q\ = Jq^ + k^-2qx. (6.27) 

The initial vacuum state is recovered when rT^ — > implying that the particle content of the 
initial state vanishes. In the latter limit, in particular, we have 

jim A(a) = 0, (6.28) 

lim B.(a) = iJdAnkef'^n,,^-'-nrk-m 

^ ' 2 |/rflnA;fc3-2"|2 V J 



^■^The result will be given in the case of the gravitons. The results for the phonons can be obtained from 
the ones of the gravitons by replacing a — /3. 
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the same result holds, as previously remarked, by exchanging a withe /3. By introducing the 
rescaling y = q/k the explicit form of Eq. fl6.29p becomes 



3 / In A; A;^"^" / In g g3-2a ji^ 



BJa) 



dx 



2 |/rflnA;A;3-2"r 
After integration over the x variable, the numerator of Eq. fl6.3Up becomes 

dx 



(6.30) 



dlnkk 
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3-2a 



dhiqq 
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1 (1 + ?/2 _ 2yx)°' 



dink r 
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dy y 



1-2q 



2(1-q) 



\i-y 



2(l-a) 



(6.31) 



The value of the integral over y = q/k can be estimated as 



4(1 



1 



3 - 2a 



1 - 

ymin 



+ 



3 -4a 



y3-4a _ ^ 
ymax 



(6.32) 



The limit T/max — ^ oo can be takeio since (3 — 4a) < (and (3 — 4/3) < 0). Conversely the 
value of T/min = qo/k forbids taking the limit j/mm — )■ since (3 — 4a) and (3 — 4/3) are both 
negative as established in Eq. fl6.15p . Since the term containing ?/min dominates the integral 
of Eq. (I6.3ip can be estimated up to sub-leading corrections. Integrating over k and going 
then back to Eq. fl6.24p we have 



lim g^'^\r,T) 

|fcr|->0, |fcr|->0 



3. 



(6.33) 



The value given by Eq. fl6.33p coincides with the value obtained, in the case of a single 
degree of freedom (and for the squeezed vacuum state) as discussed in section [2l As men- 
tioned in section [2] the normal ordering in the operator (or its absence) does not affect the 
degree of second-order coherence as long as the average number of particles per field mode 
is much larger than one. This is, a posteriori, exactly the physical limit of Eqs. f l6.1ip - 
fl6.12p . In connection with Eq. fl6.33p two final remark are in order. As mentioned after Eq. 
fl6.22p the limits appearing in Eq. fl6.33p signify that the result holds for typical comoving 
three-momenta larger than the Hubble radius; the strict mathematical limit would instead 
stipulate that r — > 0. The limit defined in Eqs. f l6.22p and fl6.33p demands that the integrals 
appearing in the numerator and in the denominator of the degree of second-order coherence 
are evaluated for typical wavelengths larger than the Hubble radius. 



^^It is worth stressing that the negative sign of the combinations (3— 4a) and (3—4/3) is a direct consequence 
of the determinations of cosmological parameters reported in Eqs. (|6.ip . (|6.2I) and (|6.3p . 
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6.3 Classical stochastic variables 



The result obtained in Eq. f l6.33p imphes that the rehc gravitons and rehc phonons are highly 
bunched and their statistics is super- Poissonian. The large-scale curvature fluctuations might 
also be classical stochastic variables. In the latter case, however, the degree of second-order 
coherence in the zero time-delay limit will typically be between 1 and 2, i.e. 1 < g^"^^ < 2. 
In this case the averages appearing in g^"^^ will have to be interpreted as stochastic averages. 
The lower limit (i.e. g^"^^ = 1) is easy to justify. Suppose, indeed, that the relic phonons (or 
gravitons) are just produced independently. This means that the intensity X can be viewed 
as a classical (discrete) variable characterized by a Poissonian probability distribution. In 
this case (X'~) = (X)'' and, consequently, g^'^^ = 1. In is also possible to conceive a completely 
classical situation where the source is chaotic in such a way that (X'') = r! (X)^. In the 
latter case the probability distribution for the (possibly time- or space-dependent) intensity 
can be written as -P(X) = X ^ exp [—X/X], where X = (X). 



6.4 Finite volume effects 

As a next step of complication it is useful to analyze the case where the number of particles 
in each Fourier mode is the same but second-order correlations are allowed in the initial 
state, i.e. 

nk = n, Ck = C- (6.34) 

Following the same steps outlined in the vacuum situation, the analog of the limit given in 
Eq. (Km is 

hm ,(%r.)-3 I (3-2a)^ 
IkriJ'fkri^o^ ^"^'^^--^^ klV C (2n + l)2 4a-3 • ^^"^^^ 

Concerning this expression few comments are in order: 

• the first term (i.e. 3) remains also in the limit n (i.e. in the vacuum case); 

• the second term goes to zero in the infinite volume limit (which is the one which must 
be enforced) and also in the case of a Bose- Einstein distribution (i.e. C 1); 

• the numerical factors in the second term depend upon the slow-roll parameter. 



The result of Eq. fl6.35p shows, in a specific example, that second-order correlation effects 
possibly present in the initial state vanish when the volume goes to infinity and the average 
multiplicity goes to infinity (while their ratio is kept fixed). The implications of this result 
are interesting per se but a closer scrutiny goes beyond the aims of this script. 



31 



6.5 Bose-Einstein occupation number 

As suggested by the general equations derived in section [6l the degree of second-order co- 
herence is sensitive to the overall duration of the inflationary phase even in the case Cfc = 1 
where, by definition, volume effects are absent and the generalized statistical ensemble of 
Eq. fl5.9l) reduces to the thermal (or chaotic) one (see Eq. fl4.8p ). Resorting to the quantum 
mechanical analogy discussed in section [2] it would be tempting to conclude on the basis of 
Eq. (12.111) . that the degree of second-order coherence should still equal 3. Indeed, the limit 
~ n ^ 1 and C ~^ 1 implies, from Eq. (12. lip , that g^"^^ — ?■ 3. This conclusion is incorrect 
as long as, in the realistic field theoretical case, Nk 7^ fik] more specifically 

= e^AT-l ' = ^' ^^-^^^ 

where is the (comoving) thermal momentum whose explicit value can be usefully expressed 
in Hubble units (i.e. units of the Hubble rate Hq) as 

i/o" UA2.43XIO-9J U.15xl0-9j UoA ^ ^ 



In Eq. (I6.37P the comoving value of the thermal momentuno depends upon the ratio {T/H) 
which measures the ratio between the temperature and the Hubble rate when Atot — A^max, 

'T\ _ /106.75\ 1/4/2.43 X 10-9\i/V0.01\ 

where Qp denotes the total number of spin degrees of freedom for T > 200 GeV. If Atot ^ 
Amax, will bccome arguably much smaller than Hq as Eqs. (I6.37p and (I6.38p imply 
immediately. 

The uncertainty in the estimates of Eqs. (I6.37P and (I6.38P resides not only in the (in- 
evitably) unknown value of the total number of e-folds but also in A^max whose value cannot 
be precisely assessed even within the consistent lore provided by the conventional inflation- 
ary scenarios. Indeed, the uncertainty affecting the determination of Amax is due to the lack 
of a specific knowledge of the post-inflationary thermal history. Suppose, for instance, that 
right at the end of the inflationary phase, the standard radiation-dominated phase starts. 
In this case we have that 

ATm,. = 62.2 + ^\ - In f - In f \ (6.39) 

2 VlO-V V0.7/ 4 V4.15xl0-5; ^ ' 

Recalling that, according to the WMAP 7 yr data, the amplitude of the scalar power spec- 
trum at the pivot scale fcp = 0.002 Mpc~^ is given by A-r = 2.43 x 10^^, the estimate of Eq. 
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As already mentioned we shall normalize to 1 the present value of the scale factor qq = 1. 
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fl6.34p becomes (see, e.g. [55] ) 

iVmax = 63.6 + ^lne. (6.40) 

The figures given in Eqs. (16.341) are just a dim indication since A^max can indeed be much 
larger. Suppose, for instance, that right after inflation the Universe expands at a rate which 
is slower than radiation. In this case A^max increases. In particular, if, after inflation, the 
energy density of the plasma is dominated by a stiff source with sound speed coinciding with 
the speed of light we get to the estimate 

A^max = 78.3 + ^lne, (6.41) 

where it has been assumed that the stiff phase starts right after inflation and stops right be- 
fore big-bang nucleosynthesis (see, e.g. [23 EH ESI ES] and references therein). By definition 
A^max is derived by requiring that the present size of the Hubble radius is all contained in 
the event horizon at the onset of the inflationary phase. Being optimistic we can say that 
A'max = 63 ± 15 which is anyway a pretty large indetermination. The indetermination on 
the specific values of A^tot and A^max implies that, unless A^tot — A^max, the thermal wave- 
length k^^ will be much larger than the present value of the Hubble rate. Bearing in mind 
the previous caveats, using Eq. (I6.36P and recalling Eqs. (I6.24p -( l6.26p we shall have that 
B,{a) = \^{a)\y\V{a)\^ where 

ma) ^ v/3/..,.««-/d,,'-coth (^) |;";.'-eoth {^) (6,42) 

V{a) = V2j d\nk coth . (6.43) 

The evaluation of the integrals can be performed with different methods and even numeri- 
cally. It is instructive, however, to derive an explicit analytic estimate based on the observa- 
tion that cothx can be approximated with 1/x for x < 1 and with 1 for x > 1. The integrals 
of Eqs. (I6.42P and (I6.43P can therefore be evaluated by using this simple approximation 
scheme which can be improved by keeping further terms in the expansion. To lowest order 
we will then have that 



Via) 



1-a) 



3 - 2a 



^a^ - (2^t)^-^" 



(6.44) 



The leading and subleading terms in Eq. (I6.44p are determined by the hierarchy between 
/cmin, Hq and /ct- In particular, if fcmin > 2A;t then it will always be true that k > 2kT 
since, by definition of /cmim k cannot be smaller than fcmin- But then from Eq. fl6.36p and 
(I6.43P (2nfc + 1)^1 for all the range of the momenta. This situation happens, in particular, 
when A^tot ^ A^max- In this situation the degree of second-order coherence will reproduce 
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the vacuum case, i.e., for \kr\ ^ 1 and |A;r| ^ 1, g^'^\r,T) — )■ 3. In the opposite situation 
kmin > 2/ct but then the leading term in Eq. f l6.44p will be the one coming from krain (which 
will be assumed to coincide with Hq for the purpose of numerical estimates). The rationale 
for the latter statement stems from the value of (3 — 2a) which is always negative in the 
case of conventional slow-roll dynamics. The same kind of considerations can be used for 
the estimate of N'{a) which can also be written, after integration over z, as 



M{a) = Jdk fc^-^" J dyy^-^'^ coth 



( k W (2A;tA)2(i-'^) 



V2fcT/ 1(1 -2a)(2-2a) 



2(1 -a) \ k ) {l-2a) ] ^ ' 

Using Eq. fl6.14p expressing the relation between the indices a (and /3) and the tensor 
(and scalar) spectral indices, the final result for the degree of second-order coherence can be 
written, in the limit /cr — )■ and A;r — > 



^ (nt-2)(nt + 2)(nt + 8) 
2rt 

^ _ (3ns - Ant - 11) (2n, - Ut - A){nt + 2) 
8(ns - l)(^t + 3) 



(6.46) 
(6.47) 



where Eq. f l6.46p holds for the gravitons while Eq. fl6.47p holds in the case of the scalar 
phonons. Both expressions have been obtained by using the relations previously derived in 
Eq. (16. 6p and by demanding, as previously explained, that A^tot — A^max- It is interesting to 
notice that as long as and Us — 1 both diminish the second term in the square brackets 
increases above 1 and it can happen that g^'^\r,T) < 3 possibly becoming even negative. 

Let us now assume, for a moment, that the intensity correlations are experimentally 
accessible and that g^"^^ (i.e. the degree of second-order coherence in the large-scale limit) 
can be directly measured. The considerations developed in the present section suggest three 
considerations: 

• if g^"^^ = 3 the HBT correlations would imply that the fluctuations are bunched, super- 
Poissonian and probably coming from a vacuum initial state; 

• if g^"^^ > 3 the initial state contained second-order correlations which are not char- 
acterized by Bose-Einstein statistics but which, nonetheless, dominate at large scale; 
this would probably be a remnant of the initial state and would imply, within the 
inflationary lore, iVtot ~ A^max ~ N^m] 

• if 2 < g^'^'^ < 3 the initial conditions are dominated by a thermal ensemble and A^tot ~ 

' max 1 
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• if 1 < g^'^^ < 2 the intensity correlations are still super Poissonian but not squeezed. 

The possibilities listed above are the result of the preliminary discussion reported here and 
must be sharpened further. At the same time they illustrate how a direct study of the 
intensity correlations would make the problem of the large-scale initial conditions much less 
elusive. Finally, recalling the terminology introduced at the end of section [2], curvature 
phonons can even be super-chaotic. 

The value of 'g^'^^ measures, in the context of HBT interferometry, the statistical ten- 
dency of the phonons to distribute themselves in bunches rather than randomly and thus 
obeying a Poissonian distribution. If the degree of second-order coherence is larger than 1 
the positive correlation between the particles arriving at the HBT detectors is dubbed as 
super-Poissonian. In the case of the curvature phonons the degree of second-order coherence 
can even be larger than 2. In this case it is natural to talk about the possibility that the 
statistics is super-chaotic since the case g^'^'> = 2 corresponds to the statistics of chaotic (i.e. 
white) light. The bounds summarized in the previous paragraph are illustrated in Fig. [3] 





A : generic remnants 

of initial state 
B: thermal initial 

state 
C: generically 

super-Poissonian 



Super-chaotic limit 



Chaotic limit 



Poissonian boundary 



Figure 3: The degrees of first-order and second-order coherence are reported in the case of 
the scalar phonons. 

where, in the left and right plots, the degrees of first-order and second-order coherence are 
reported. The regions A, B and C correspond, respectively, to a super-chaotic state, to a 
squeezed thermal state and to a super-Poissonian statistical ensemble not necessarily related 
to quantum mechanics. 



7 HBT temperature and polarization correlations 

Some preliminary considerations will now be developed with the purpose of suggesting that 
HBT correlations could be directly assessed by studying higher order temperature and po- 
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larization correlations in the limit of large angular scales. Consider, to begin with, the 
brightness perturbations which fully describe the temperature and polarization anisotropies 
in the ACDM model neglecting, for simplicity, the tensor mode^ which are indeed absent 
in the vanilla ACDM: 

Ai(n,r) = -^^^ J d^kAi{k,fi,T), Ap(n,r) = J (i=^A; Ap(/c, /x, r); (7.1) 

as usual, fi = k-n denotes the projection of the Fourier mode on the direction of propagation 
of the CMB photo,§ nofce. furthermore, that A p(?7,,r) = AQ(?2,r). Since A±(n,r) = 
AQ(n, r)±z Au (hr), it transforms as a spin ±2 for rotations around a plane orthogonal to the 
direction of propagation of the radiation. The three-dimensional rotations and the rotations 
on the tangent plane of the sphere at a given point combine to give a 0(4) symmetry group 
[67] . Generalized ladder operators raising (or lowering) the spin weight of a given function 
can then be defined as [67l [68] : 



sm v 



(sin^9)^^ n = {^,ip). (7.2) 



In real space the E-mode and the B-mode polarization will have spin weight s = 

A^{h,T) = -^{K^}\h)[K^'\h)A^ih,T)] + Ki-'\h)[^^^^^^ (7.3) 

AB(n,r) = ^{ir«(n)[iri')(n)A+(n,r)] -iri-^)(n)[iri-')(n)A_(n,r)]}. (7.4) 

and will therefore be invariant under rotations around n exactly as Ai(n, r). Recalling that 
fi = cos?? the derivatives with respect to can be traded for derivatives with respect to /i 
and Eqs. (17. 3p and (17.41) imply, in the ACDM framework and in the absence of gravitons, 

AE(n, r) = -dliil - fi')Ap{n, r)}, AB(n, r) = 0, (7.5) 

where the second equality is easily deduced since the derivatives with respect to vanish 
in the absence of the tensor contribution of the gravitons to the brightness perturbations. 
Therefore, in the ACDM model we have at our disposal three complementary degrees of 
second-order coherence which can be defined 

(2).. . X (Ai(m, r) Ai(m, r) Ai(n, r) Ai{h, r)) 

' (Ai(m,r)Ai(m,r))(Ai(n,r) Ai(ri,r))' ^ ' 



^^The same considerations developed in this section can be however extended to the case of the tensors 
by considering higher-order correlations of the B-niode polarization, if and when they will be measured. 

^®The angular variables ^ — k ■ n defined and used in the present section have no relation with ^{x,t) 
used and defined in section [3] at Eq. p.Sp ; similarly the angular variable of Eq. ()7.2|) is not related to the 
inflaton background field. 
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(2) / - - ^ 



(2) / ^ ^ \ 



(AE(m, r) AE(m, r) AE(n, r) AE(n, r)) 
(AE(m, r) AE(m, r)) (AE(n, r) AE(n, r)) ' 
(AT(m, r) AE(m, r) AT(n, r) AE(n, r)) 



(7.7) 
(7.8) 



(AT(m, r) AE(m, r)) (AT(n, r) AE(n, r)) ' 

The degrees of second-order coherence defined in Eqs. (I7.6p . (I7.7P and (17. Sp can be connected 
with the degree of second-order coherence introduced in sections and El In the gauge defined 
by Eq. (11. ip the brightness perturbations of Eq. (17. ip obey, in Fourier space, 



drAp + {ikfi + e')Ap 



3,. 



ikficj) + e' 



Aio + fiVb + 



-Sp{k,T) 



^?)Sp{k,T) 



(7.9) 
(7.10) 



where Sp{k,T) can be expressed as the sum of the quadrupole of the intensity, of the 
monopole of the polarization and of the quadrupole of the polarization, i.e., respectively, 
Sp{k,T) = (Ai2 + Apo + Ap2); note that, in Eqs. frTOjl and fmoj) . e' and e(r, tq) denote, 
respectively, the differential optical depth and the optical depth itself 



e = XeUe a a, 



■yej 



To) 



(7.11) 



and should not be confused with e defined, in the present paper, as one of the slow-roll 
parameters (see, e.g., Eq. (16. 5p ). The line of sight solution of Eqs. (17. 9p and (I7.10p can be 
written, respectively, as 



Ai(fc,/i,ro) 



TO 



/C(r) 



Aio + </) + yU^^b + 



3/i2 



-Sp 



-ifix{T) 



+ 



Ap(fc,/i,ro 



r dre-"("'"°)(0' + V'')e"'''"^"^c^r, 
Jo 

= ^(l-/i') r )C{T)Sp{k,T)e-^'^^^-^^^dT, 
4 Jo 

£(t,to) jg ^Yie visibility function and x{t) = k{TQ - 



(7.12) 
(7.13) 



where /C(r) = s'e~^^'^''^°^ is the visibility function and x{t) = k{TQ — r). The visibility 
function can be approximated as a double Gaussian with two peaks roughly corresponding 
to the redshifts of recombination and reionization, i.e. Zj-ec — 1088.2 ± 1.2 and Zreion = 
10.5 ± 1.2 according to [HE]- The semi-analytical parametrizations of the visibility function 
(such as the ones of [69l [70]) are relevant when investigating the degree of second-order 
coherence for angular scales than the degree. In the present analysis the focus will be on 
the large-angular scales corresponding to typical multipoles i < y/z^ where the finite width 
of the visibility function is immaterial and the opacity suddenly drops at recombination. 
This implies that the visibility function presents a sharp (i.e. infinitely thin) peak at the 
recombination time. Thus /C(t) is proportional to a Dirac delta function and e"'^*^"^'^"-* is 
proportional to an Heaviside theta function. Under the latter approximations, Eq. ( ]7.9p 
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leads to the well known pair of separated contributions, i.e. the Sachs- Wolfe (SW) and the 
integrated Sachs- Wolfe (ISW) contributions: 



Ar)(A;,/.,ro)=r-^) , (7.15) 



Ai(fc, /X, To) = Ar>{k, /i, To) + Ar^'ik, /X, To), (7.14) 

% 
5 

Af'^\k, /i, To) = r (0' + ^')e-'^^'^^^ dr, (7.16) 

*^ '^rcc. 

where, by definition, x(rrec) = Vrcc- Equations (17. 15p and (I7.16P can be evaluated within 

various approximation schemes. The SW and the ISW contributions can be separately 
evaluated. In particular the ordinary SW contribution becomes 

^^\k,^L,T„) = -'Eih3ls(g,^)e--y- (7.17) 



5(9) = l + ^-^+ ^"'"":3 (7.18) 



4 16 16(v/yTT-l 
3q 3q 

while the ISW contribution is: 



Af^'^\k, /i. To) = -27^(fc, Ti) r 9.r^(r)e-^^-'W dr, (7.19) 
Tn{r) = l-'^J\\r')dr'. (7.20) 



a2(r) Jo 

Both in Eqs. (I7.17P and (17.191) . TZ{k, r) denotes the constant value of curvature perturbations 
at Ti < Tcq. By further approximating the integrand in Eq. (17.191) the whole large-scale 
contribution can be written, for the present purposes, as 



rcc 



where 



Ai(fc,ro) = n{k,n)e- 
S{qrec) = -{%^ + [5.7-^,4, (7.21) 

g,ec = — = = 3.04 f . (7.22) 

^ Oeq 2rec + l V 0.134/ ^ ' 

Equation (]7.2ip directly relates the curvature perturbations to the brightness perturbations. 
It then follows that the two point function of temperature perturbations bears the mark, 
up to time-dependent factors, of the two-point function defined in the context of first-order 
coherence effects (see section |3]). In particular we have that 

(Ai(m,ro)Ai(n,ro)) =E^^^^C'r^^(^-^), (7-23) 



Cr = Jdlnqq^ {2n, + 1) \^2N, + 1 - 2 cos ajN,{N^ + 1) 



j/[l/,ec].(7.24) 
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In a similar fashion the degree of second-order coherence can be estimated as 



(Ai(m, r) Ai(m, r) A,{h, r) Ai(n, r)) = ^(2£ + 1) 



Zi + QiPeim ■ n) 



(7.25) 



where 



Qe 



<S {qr< 



dlnq 
dlnq 



The degree of second-order coherence of Eq. (17. 6p can therefore be written as 



(2) / - - N 

(m,n,r) 



\Ei{2i + l)Cr\' 



(7.26) 
(7.27) 

(7.28) 



All the inequalities established for the degree of second-order coherence and all the consid- 
erations presented before are also applicable to Eq. (I7.28p : in the limit m ■ n — )■ 1, using the 

(2) _______ 

well known identities, gj^j^ coincides with the result expressed, for instance, by Eq. (16.231) . In 
the cases m-n ^ la specific angular dependence should be taken into account when trying to 
infer the degree of second-order coherence from the observational data. The considerations 
developed in this last section exclude the presence of the tensor modes which can be however 
included without problems. 



8 Concluding remarks 

In conventional Hanbury Brown-Twiss interferometry the statistical properties of the source 
are often part of the experimental setup but the space-time dimensions of the emitters need 
to be determined. For large-scale curvature perturbations the reverse is true: while the 
statistical properties of the source are unknown, the gross uniformity of the temperature 
fluctuations at last scattering implies that curvature perturbations prior to matter-radiation 
equality had typical wavelengths larger than the Hubble radius at the corresponding epoch. 

Can we directly scrutinize the statistical properties of the large-scale curvature perturba- 
tions without positing an excessive number of assumptions on the pre-infiationary expansion 
and on the post-inflationary thermal history? This has been the main question addressed 
in the present paper. As a partial and preliminary answer, it has been suggested that a 
useful approach to large-scale curvature perturbations of quantum mechanical origin con- 
sists in scrutinizing their large-scale coherence properties. Using then the analogy with a 
similar problem arising in quantum optics, various results have been obtained and they can 
be summarized as follows: 
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• in the limit of wavelengths larger than the Hubble radius the degree of first-order 
coherence goes always to 1 in spite of the correlation properties of the initial state; 

• the degree of second-order coherence bears the mark of the statistical properties of the 
initial state; the curvature phonons are bunched and their degree of bunching exceeds 
the typical value of a chaotic source; 

• direct limits on (or specific determinations of) the degree of second-order coherence 
from temperature and polarization maps can probe the correlation properties of large- 
scale gravitational fluctuations; 

• the degree of second-order coherence does depend, in a computable manner, upon the 
values of the slow-roll parameters, upon the nature of the initial state and upon the 
duration of the inflationary phase ; 

• a set of model-independent limits on the degree of second-order coherence has been 
derived in the form of a collection of inequalities which can be tested explicitly once 
the degree of second-order coherence is defined in terms of the correlators involving 
the relevant brightness parturbations. 

On a more technical ground, the tenets of the quantum theory of optical coherence have been 
carefully translated to the quantized treatment of the scalar and tensor normal modes of the 
geometry. The correct quantum mechanical definition of the degree of second-order coherence 
has been derived and it has been shown to be equivalent, in the limit of a large number of 
phonons per Fourier mode, to the standard normal-ordered definition customarily employed 
in Hanbury Brown- Twiss interferometry and based on the quantum theory of photoelectric 
detection. 

Pre-infiationary initial conditions are often assigned (or imposed) by combining theoret- 
ical prejudice with loose elements of phcnomcnological consistency such as the suppression 
or the increase of large-scale power spectra. Instead of arguing that the theoretical prej- 
udice necessarily selects a specific number of infiationary e-folds, a preferred set of initial 
conditions, a unique pre-infiationary history it is also worthwhile to take a more modest 
approach and to ask ourselves wether it is possible to gain informations on the correlation 
properties of pre-decoupling initial conditions by using the logic of Hanbury Brown-Twiss 
interferometry, i.e. the study of the intensity correlations of the scalar and tensor fluctu- 
ations of the geometry. The results reported in the present analysis represent a flrst step 
along this direction. 
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